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>• ! Abstract: A general approach to a fast simulation of electromagnetic showers 
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using parameterizations of the longitudinal and radial profiles in homogeneous 
O ' and sampling calorimeters is described. The dependence of the shower devel- 

^ opment on the materials used and the sampling geometry is taken into account 

. explicitly. Comparisons with detailed simulations of various calorimeters and 

O I with data from the liquid argon calorimeter of the HI experiment are made. 
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1 Introduction 

• ^ 

X. 

^ , In calorimeter simulation different tasks can be distinguished: calorimeter studies, physics 

analysis, and feasibility studies. A detailed simulation, where all secondary particles are 
tracked individually down to some minimum energy and where the response is predicted 
from "first principles", is required for accurate calorimeter studies. For physics analysis 
and feasibility studies large number of Monte Carlo events may have to be produced. Us- 
ing individual particle tracking, the computing time needed for such kind of simulations 
increases approximately linear with the energy absorbed in the detector and can easily 
become prohibitive. Using parameterizations for electromagnetic (sub) showers can speed 
up the simulations considerably, without sacrificing precision. The high particle multiplic- 
ity in electromagnetic showers as well as their compactness and the good understanding 
of the underlying physics makes their parameterization advantageous. 

Using an Ansatz by Longo and Sestili a simple algorithm for the description of 
longitudinal shower profiles has been used successfully for the simulation of the UAl 
calorimeter [El. Later, this Ansatz has been extended to the simulation of individual 
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showers, taking their shower-to-shower fluctuation and correlations consistently into ac- 
count P, H- For the parameterized simulation of radial energy profiles no conclusive 
procedure has been established until now. 

In homogeneous media, a scaling of the longitudinal and radial profiles in radiation 
lengths and Moliere radii respectively does not lead to a material independent description 
of electromagnetic shower development. In sampling calorimeters, the shower shapes 
depend in addition on the sampling structure. We have extended the above Ansatz for 
parameterized simulation of longitudinal profiles by taking the material and geometry 
dependence of the parameters into account and developed a new algorithm to simulate 
radial energy distributions [§. Correlations between the longitudinal and radial shower 
development have been included. 

2 Procedure 

To arrive at a general description of electromagnetic shower development, we performed 
detailed Monte Carlo simulations, on a grid of 1.0 Xq in depth and 0.2 Moliere radii 
laterally, for various homogeneous media and sampling calorimeters, using the GEANT 
package 0]. The materials used were Cu, Fe, W, Pb, U, and scintillator and liquid 
argon. In a first step only average shower profiles in homogeneous media were analyzed, 
from which scaling laws for the material and energy dependence of the parameters have 
been extracted. Starting from the relations which describe the average behavior of the 
parameterized quantities, we developed parameterizations for individual electromagnetic 
showers in homogeneous calorimeters, taking fluctuations and correlations into account. 

The parameterizations in homogeneous media are a first approximation for electromag- 
netic shower development in sampling calorimeters, which are viewed as consisting of one 
single effective medium. The inhomogeneous material distribution in sampling calorime- 
ters influences however the exact behavior of the shower shapes. This is explained mainly 
by the transition effect which depends on the shower depth H, These effects have 
been taken into account by adding geometry dependent terms to the parameterizations 
for homogeneous media, which can be easily calculated from the sampling geometry. 

3 Parameterization Ansatz 

The spatial energy distribution of electromagnetic showers is given by three probability 
density functions (pdf). 



describing the longitudinal, radial, and azimuthal energy distributions. Here t denotes 
the longitudinal shower depth in units of radiation length, r measures the radial distance 
from the shower axis in Moliere units, and is the azimuthal angle. The start of the 
shower is defined by the space point, where the first electron or positron bremsstrahlung 
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process occurs. A gamma distribution is used for the parameterization of the longitudinal 
shower profile, f(t). The radial distribution, f(r), is described by a two-component Ansatz. 
In 0, it is assumed that the energy is distributed uniformly: = l/27r. 



3.1 Longitudinal shower profiles — homogeneous media 

It is well known that average longitudinal shower profiles can be described by a gamma 
distribution [|I|: 

lldE{t)\ _ _ m-^PeM-Pt) ..^ 
\E dt I ~ ' ~ V{a) ■ ^ ' 

The center of gravity, (t), and the depth of the maximum, T, can be calculated from the 
shape parameter a and the scaling parameter (3 according to 

W = ^ (3) 
T = — . (4) 

Longitudinal electromagnetic shower development in homogeneous media had been stud- 
ied analytically by Rossi [jlOl- An important result of the calculations using "Rossi Ap- 



proximation B" is that longitudinal shower moments are equal in different materials, 
provided one measures all lengths in units of radiation length (Xq) and energies in units 
of the critical energy {Ec). Numerically, Ec can be calculated according to [|Tl|] 



Z 



1.1 



Ee = 2.66(Xo-) . (5) 



For the depth of the shower maximum 



T oc Iny = In — (6) 

Er 



is predicted Jiy . 



It is therefore desirable to use T in the parameterization. This is demonstrated in 
Fig.|l], where the average depth of the shower maximum for various homogeneous mediaQ, 
Thorns is plotted versus in the energy range from 1 to 100 GeV. As a second variable 
a is used. In this case the parameterization depends on the charge number Z of the 
medium, as can be seen in Fig.^. The lines in both figures correspond to fits to GEANT 
simulations according to 

Thom = Iny + ti (7) 
ahom = fli + (a2 + cts/^) Iny. (8) 



^ The index "horn" in the following formulae indicates the validity for homogeneous media. For 
sampling calorimeters the index "sam" will be used. 
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The values of the coefficients are given in Appendix, where all formulae and numbers, 
which will be given in the following, are summarized. 

Assuming that also individual profiles can be approximated by a gamma distribution, 
the fluctuations and correlations can be taken into account consistently (for details refer 
to P). For each single GEANT-simulated shower, T and a are determined by fitting 
a gamma distribution. The logarithms of T and a are used for the parameterization 
since they are found to be approximately normal distributed. For the parameterization of 
(In T/iom) and (^v^ahom) the logarithms of equations ^ and ^ are used. The ?/-dependence 
of the fluctuations can be described by 

o = {si + S2lny)-\ (9) 

The correlation between In Thom and In ahom is given by 

p{\nThomA'^ahom) = P = ri + rslny. (10) 

The dependence of these quantities on y is shown in Fig.|3 for various materials together 



with the parameterizations (see Appendix A. 1.2 ). 



From these formulae, correlated and varying parameters a, and Pi are generated ac- 
cording to 



with 

C 



a(lnT) 
a(\na 




and Pi = (aj — 1)/Tj and zi and Z2 are standard normal distributed random numbers. The 
longitudinal energy distribution is evaluated^ by integration in steps of At = tj — = 

1^0, 

Jtj-i ^ycni) 

It is worthwhile to mention that only one of the five quantities needed, {In ahom), depends 
explicitly on the material, while for the other four this dependence is absorbed by using 
y instead of E. 

In Fig.^ longitudinal proffies of GEANT and parameterized simulations for a lead 
glass calorimeter (SF5) are compared. Shown are the mean proffies and the mean + 1 
RMS in each Xq interval. While the means are in perfect agreement, the fluctuations 
are underestimated by the parameterized simulations at low energies, indicating that the 
description of individual profiles by gamma distributions becomes a worse approximation 
with decreasing shower energy. Comparisons for other materials (Fe, Cu, W, Pb, U) 
are of comparable quality as those in Fig.^ 0. In the next sections we will show, how 
the sampling fluctuations in sampling calorimeters can be used to improve the shape 
fluctuations at low energies. 

^ The GAMDIS function of the CERN computer hbrary is used. 
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3.2 Sampling fluctuations 



In fast simulations, sampling calorimeters consisting of a complicated but repetitive sam- 
pling structure are usually described by one single effective medium (the formulae to 
compute effective material parameters are summarized in Appendix [A.2.1|) . The sampling 
fluctuations, the scaling of the deposited energy to the visible energy using an appropri- 
ate sampling fraction, and the effects of the sampling structure have to be considered in 
parameterized simulations explicitly. 

The simulation of sampling fluctuations are done conveniently with a gamma distri- 
bution: 

G(a,6) = ^-f^ (12) 
r(a) 

with 

(x) = ^ , a\x) = ^. (13) 
The energy in each longitudinal integration step, dE{t), is fluctuated^ according to equa- 



tion |T2| choosing 

a = — ^ and b = —. (14) 

c c 

It is then easy to show that the central limit theorem will ensure the total energy to be 
normal distributed obeying the usual formula for the sampling fluctuations: 

Using this procedure the occurrence of negative energies is automatically avoided. Addi- 
tional fluctuations of the longitudinal shape are introduced, leading to a better agreement 
in the shape fluctuations. This method is also used to fluctuate energy depositions of real 
particles (electrons, hadrons), when they are tracked individually through an effective 
homogeneous volume. 



3.3 Longitudinal shower proflles — sampling calorimeters 

The inhomogeneous material distribution in sampling calorimeters influences the exact 
behavior of the shower shapes. In the flrst stages of electromagnetic shower development 
the signal is dominated by electrons and positrons. Behind the shower maximum low en- 
ergetic photons become more and more important. The transition effect, being explained 
mainly by the absorption properties of low energetic photons, must in turn depend on 
the shower depth. Consequently, the signal ratio of electrons to minimum ionizing par- 
ticles, e/mip, decreases continuously as the shower propagates longitudinally. Thus the 
signal maximum in a sampling calorimeter occurs at an earlier depth than expected for 
a homogeneous calorimeter with the same effective material properties. This can be seen 
from Fig.^ (left upper corner), where (InT) for homogeneous media is compared to the 

^ The RANGAM function of the CERN computer hbrary is used. 
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values in five different sampling calorimeters. In addition, the amount of the shift of (In T) 
depends on the exact geometrical arrangement. 



The parameterization of the longitudinal shape as given in section |3.1| for homogeneous 
media can therefore not be used for sampling calorimeters directly. Instead it may be 
understood as a first approximation to which geometry dependent corrections have to be 
added. We use the sampling frequency 



(16) 



and the value of e/mip (averaged over the shower depth) to account for the shower depth 
dependence of the transition effect, da and dp denote the thickness of the active and 
passive layers, respectively. If e/mip is not known, a sufficiently good approximation for 
many calorimeters [|1^] with charge numbers Zp and Za is given by 

(17) 



1 + 0.007(Zp - Za) mip 

Averaged over the whole shower, e/mip remains energy independent for E^l GeV. 
The average longitudinal profiles can now be parameterized according to 

Tsam = Thorn + tlFg ^ + t2{l — e) (18) 
"sam = CHhom + aiFg^, (19) 

and the quantities used for the simulation of individual showers are given by 

(InT,,^) = \n{exp{{lnThom))+tiFs^ + t2{l-e)) (20) 
{\nasam) = In (^exp((lna,,o^)) + aiF^^) . (21) 

The fluctuations, a{\nTsam), o'(lnasam) and the correlation, p{\ia Tsam sam), are de- 
scribed with the help of the same formulae as in the case of homogeneous media (see 
Appendix and [AlTSl) . 



Fig.|^ summarizes the parameterization for sampling calorimeters. The expectation 
value of InT no longer scales with y. The expectation value of In a depends on the 
material and the sampling geometry. The fluctuations and correlations of the parameters 
can still be approximated without any explicit material or geometry dependence. 

In Figs.|^ to 1^ GEANT and parameterized simulations of the lead liquid argon calorime- 
ter (IFE) of the HI experiment [0, |^ are compared. The GEANT simulations were per- 
formed with low energy cuts (e-cut= 200 keV, 7-cut= 10 keV) and a detailed geometry 
description, including for example copper pads and GIO layers. These simulations were 
not used to tune the parameterizations. Both, average longitudinal profiles and their fluc- 
tuations (including sampling fluctuations) are in very good agreement (see Fig.|^). The 
energy containment (see Fig.^ and the energy resolution (see Fig.|D as a function of the 
longitudinal calorimeter length are also well predicted. Comparisons with detailed simu- 
lations of other calorimeters (Fe-LAr, Cu-Sc, W-LAr, Pb-LAr, U-Sc) show a comparably 
good performance [Q. 
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3.4 Radial shower profiles — homogeneous media 

Average radial energy profiles, 

1 dE(t,r) , , 

at different shower depths in pure uranium are presented in Fig.^ These profiles show a 
distinct maximum in the core of the shower which vanishes with increasing shower depth. 
In the tail (r ^^li?^,/) the distribution looks nearly fiat at the beginning (1 — 2Xo), becomes 
steeper at moderate depths (5 — GXq, 13 — 14Xo), and becomes fiat again (22 — 23Xo). 
A variety of different functions can be found in the literature to describe radial profiles 
1^ , |15|, |16|, 0, |1^, 1^. We use the following two component Ansatz, an extension of 0: 



fir) = pfc{r) + {l-p)fT{r) (23) 
2ri?^ , , 2rRl 



with 

< p < 1. 

Here Rc (Rt) is the median of the core (tail) component and p is a probability giving the 
relative weight of the core component. For the shower depth 1 — 2Xq the distributions 
/(r), pfc{r), and (1 — p)fT{r) are also indicated in Fig.||. 



The evolution of Rc, Rt, and p with increasing shower depth is shown in Fig.[Ty for 
100 GeV showers in iron and uranium. We use the variable r = t/T, which measures 
the shower depth in units of the depth of the shower maximum, to generalize the radial 
profiles. This makes the parameterization more convenient and separates the energy 
and material dependence of various parameters. The median of the core distribution, 
Rc, increases linearly with r. The weight of the core, p, is maximal around the shower 
maximum, and the width of the tail, Rt, is minimal at r ^ 1. This behavior can be 
traced back to the radial profiles shown in Fig.^. 

The following formulae are used to parameterize the radial energy density distribution 
for a given energy and material: 

Rc,hom{T) = Z1+Z2T (24) 

RTMmiT) = ki{exp{k3{T - k2)) + exp(fc4(r - ^2))} (25) 
Phom{r) = pi exp <^ — — - - exp I — — - I \ (26) 

[ P3 \ P3 J } 

The parameters zi - ■ - p^ are either constant or simple functions of InE or Z (see Appendix 
|A.1.3| for details). The complicated evolution of Rt and p with the shower depth and the 



dependence on the material can be explained mainly with the propagation of low energetic 
photons 0. The offset in Rj- between iron and uranium (Fig.|10D for example, indicating 
a wider distribution in iron, reflects the difference in the mean free path, which for 1 MeV 
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photons is approximately twice as long in iron as in uranium, if lengths are measured in 
Moliere units. 

We found a good agreement of mean radial profiles between parameterized and detailed 
simulations in Fe, Cu, W, Pb, and U absorbers for energies between 0.4 and 400 GeV. This 



is demonstrated in Fig.|T2|, where radial profiles in various shower depths are compared 
for 40 GeV showers in lead and 100 GeV showers in uranium. 

The introduction of radial shape fluctuations has to be considered with some care. 
Even if no fluctuations of /(r) are simulated explicitly, the radial energy profile at a given 
shower depth will fluctuate, because the shower maximum T and thus r varies from shower 
to shower. Another source of radial fluctuations arises from the method, which we have 
adopted for the simulation of radial distributions. The energy content of a longitudinal 
interval of length 1 Xq, dE{t), is calculated from the actual longitudinal energy density 



distribution as described in section ^A\ . This energy is divided into Ns{t) discrete spots 
of energy Es = dE(t)/Ns(t), which are distributed radially according to /(r) using a 
Monte Carlo method. This can be done easily since the pdfs, fc{f) and fri^), can be 
integrated and inverted: 

F-Hu) = uf^. (28) 

Random radii are generated according to /(r) in the following way, using two normal 
distributed random numbers Vi and Wi: 




This method leads to additional fluctuations in the energy content of every radial 
interval which follow a binomial distribution. Thus, the relation 



describes the contribution to radial shape fluctuations produced by the Monte Carlo 
method in each longitudinal integration interval. Here e denotes the energy in a given 
radial interval at a given shower depth: 

/ \ - r ff w dE{t,r) 

We investigated the possibility to tune Ns{t) in each longitudinal interval to match the 
radial shape fluctuations observed in detailed GEANT simulations^. As an example, the 

^De Angelis et al. have used a similar method to reproduce shape fluctuations |ig[| . 
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quantity cr^(e)/((e) (1 — (e))) at t = 5 — 6X0 is displayed in Fig.|TT| for detailed simulations 



and parameterized ones without any radial shape fluctuations. The difference of these 



curves, which is also shown in Fig.|TT], is approximately constant and determines A^^ in 
equation ^ (note that the variance is additive). We found that a constant contribution 
to cr^(e)/((e)(l — (e))) can be used to match the total radial shape fluctuations to a good 
approximation at all shower depths. 

Summing Ns(t) over all shower depth, the total number of spots, Nspot, needed for 
one shower can be obtained and parameterized according to 

Ns,ot = 93HZ)E'-'''. (31) 

To flnd the number of spots for each longitudinal integration interval, the density distribu- 
tion 1 / Nspot dNs{t) / dt in Fig.^is parameterized. It is described by a gamma distribution 
with parameters, which are given by the corresponding longitudinal energy proflle: 

Tspot = Tw(0.698 + 0.00212Z) and (32) 
aspot = aw(0.639 + 0.00334Z). (33) 

The total fluctuations obtained with this method are compared in Fig.|12| by adding 1 
RMS to the mean proflles. 

Additional correlations between longitudinal and radial shower development are taken 
into account by introducing a correlation between the radial pdfs and the actual center 
of gravity. 



of an individual shower. This is done by replacing r in equations 24, 25, and by Xj: 



^ t_ ^ t exp((lna)) 

T ^ ' (t)iexp((lna))-r ^ ' 



The need to introduce these correlations is demonstrated in Fig.|T3|, where integrated 
radial proflles are shown, which were calculated by summing over all longitudinal layers. 
Note that the mean integrated proflles, 

/ 1 dE{r) \ 



E dr 



are independent of energy, which is well reproduced by the parameterized simulation. The 
relative fluctuations of these distributions. 



X ctrms 

^^^^ = ^nm 

E dr 



are shown using both, r and Xj, in calculating the radial proflles. Only the simulations 
using Ti are able to predict the fluctuations observed with GEANT correctly. 
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For clarity, we summarize the steps of our algorithm as follows: Determine the energy 



dE{t) within one longitudinal integration interval as described in section |3.1| . In case 
of sampling calorimeters apply sampling fluctuations on dE{t). Evaluate the number of 
spots needed to reproduce radial shape fluctuations in this interval according to 



tj-i 



Distribute the spots with energy Es = dE(t)/Ns(t) radially according to /(r) as described 
above and uniformly in (p and in the longitudinal interval At. Finally transform the spot 
coordinates {Es,t[Xo],r[RM], (p) into the detector reference system {Es,x,y, z). 

3.5 Radial shower profiles — sampling calorimeters 

The influence of the exact geometry on radial energy profiles is rather small. At the 
start of the shower the profiles look a bit smoother than in homogeneous media. With 
increasing shower depth they approach the shapes that are expected for homogeneous 
media with the appropriate effective material. These small deviations have been taken 
into account by the following corrections to the mean profiles: 

Rc,sam = Rc,hom + - e) + Z2Fs^ exp{-Ti) (35) 

RT,sam = RrMm + ki{l - c) + k2Fg^ exp(-ri) (36) 

Psam = Phom + (1 - e)(pi +P2^S^exp(-(ri - 1)^)) (37) 



using again the sampling frequency Fs and e/mip (see Appendix |A.2.4|) . 



The total number of spots needed to simulate the radial shape fluctuations is much 
smaller than in the case of homogeneous media and no longer depends sensitively on the 
materials used. Instead, the spot number can be parameterized by 

where c measures the sampling fluctuations according to 

a c 



E ^/E 

(see Fig.0). The density distribution of the spot numbers is given in analogy to the 
homogeneous media by: 

Tspot = T,„„(0.831 + 0.0019Z) and (39) 
aspot = «.am(0.844 + 0.0026Z) (40) 

GEANT and parameterized simulations of mean radial profiles and their relative fluc- 
tuations, 

dE{t) dr 
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are compared in Fig.|T5| and Fig.|l^ for the HI liquid argon calorimeter (IFE) for vari- 
ous energies. The influence of radial leakage on containment and energy resolution is 
demonstrated in Fig.17 and Fig.[l^. The energy independence of the energy contained in 
a cylinder of radius r is well reproduced by the parameterized simulations. The energy 
resolution as defined in Fig.|T8] does not depend on radial leakage. As can be seen, this 
is correctly predicted by the parameterized simulation, when the correlation between the 
longitudinal and radial shower development is taken into account (by using Tj). 



4 Comparison with data 

We have compared parameterized simulations with test beam data from the HI calorime- 
ter, which is made of lead and liquid argon in the electromagnetic sections |1T3|, |6[|. Modules 
of the inner forward (IFE), the forward barrel (FBI), and the central barrel (CB2/CB3) 
calorimeters have been studied. Electron beams in the energy range between 5 and 80 
GeV entered the stacks under angles of 11° in the IFE and CBS, and under 35° in the 
FBI calorimeter in a test set-up at CERN. 

The energy resolution of the data can be described by 



(E) \ (E) ^ 




(42) 



Here c refers to the sampling fluctuations, b considers the noise, and cr{p)/p denotes the 
momentum resolution of the beam. In the Monte Carlo the momentum resolution was 
simulated explicitly. The electronic noise was taken into account by adding random trigger 
events to the simulated cell energies. The constant c, which is approximately 11% for all 
modules, was used to simulate the sampling fluctuations. 

The simulations were carried out with the HI detector simulation program HIFAST 
pO| , |2l|| . The algorithms described so far are part of this program, which is used for the 
mass production of Monte Carlo events in the HI detector at the HERA collider at DESY. 
To keep the required high precision of the parameterization also in complicated detector 
regions (cracks for example), the following has to be considered. If a shower develops 
partly inside cracks between adjacent modules, which in general cannot be approximated 
by a single effective medium, parameterizations will in general fail to reproduce mea- 
sured signals. In H1FAST0 we therefore do not parameterize showers, if they cross such 
boundaries. Only electromagnetic showers and sub-showers from hadronic interactions 
are parameterized which fit into one single stack. 

During analysis, a 3a noise cut was applied to both the experimental data and HIFAST 
data at the cell level, and energy clusters were built from cells containing energies above 
threshold. Energy distributions of the clusters with maximum channel numbers are com- 



pared in Fig.lT9| for all three modules considered. In addition, the energy in all other cells. 



not belonging to the selected clusters, are also shown. 

stand alone version (called GFLASH 1.4) running with GEANT and covering the same functionality 
is available for distribution. Please contact one of the authors. 
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Longitudinal profiles are shown in Fig]2^ for various energies in the IFE calorimeter. 
The mean profiles as well as the fluctuations are nearly indistinguishable between data 
and HIFAST. Energy distributions in individual longitudinal layers of CBS are compared 
in Fig.^ for 30 GeV incident electron energy, showing that not only the means and fluctu- 
ations but also the shape of the distributions are predicted correctly by the parameterized 
simulations. 



Fig.|2^ compares lateral profiles in different shower depths in the IFE calorimeter at one 
energy, and in Fig.|2^ lateral profiles in the FBI calorimeter, summed over all longitudinal 
sections, are shown for various energies. There is good agreement in the peak distributions. 
The tails of the profiles are dominated by electronic noise. 

As shown so far, parameterized simulations can predict measured calorimeter signals 
very precisely, if the shower development is confined within one single calorimeter stack. 
Using the concept of partial parameterization as described above, the influence of cracks on 
the measured signal can be reproduced as shown in Fig.|4|. We have used test beam data 
scanning the crack between CB2 and CBS, which consist of two electromagnetic (CB2E, 
CB3E) and two hadronic stacks (CB2H, CB3H). The width of the crack is approximately 
1 cm. Shown are the energies in the electromagnetic modules {Ecb2E, Ecbse), the sum 
of both {EcBE = EcB2E + Ecb3e), and the sum measured in the electromagnetic and 
hadronic modules {Eqb = Ecbe + Ecbh) as a function of the beam impact position. All 
energies are normalized to -E+20 = EcBi^caio = 20cm). The energy lost while scanning 
the crack with a 30 GeV test beam extends to about 40%, if only the electromagnetic 
sections are considered, and is still around 20% if the hadronic modules are added. The 
agreement between data and partial parameterization is quite satisfactory. 

Of the various comparisons which were made P, only a limited number is presented 
here. Other properties, which are relevant for physics analysis with the HI detector, like 
e/vr separation, were studied R |21| and confirm the applicability of the fast simulation. 



5 Timing 

The CPU time reduction depends on the complexity of the geometry description and the 
cut off parameters in the detailed simulation as well as on the type of simulated event. 
Fully parameterized simulations of electromagnetic showers in a simple (box) geometry are 
about 7000 times faster at 100 GeV (900 at 1 GeV) compared with GEANT simulations 
of a detailed geometry and with low energy cuts (e-cut= 200 keV, 7-cut= 10 keV). 

In the framework of the HI simulation program, partial parameterization of electro- 
magnetic showers is performed as described above, together with individual tracking of 



hadrons and termination of low energy particles (see also [20, 0). The gain factors for 
30 GeV showers in the HI detector (including detailed simulations of tracker volumes) 
are 200 for electrons and 25 in case of hadronic showers. Medium energy cuts (e-cut= 1 
MeV, 7-cut= 200 keV) were used in the corresponding detailed simulations. Complete 
detector simulations of HERA events (ep scattering at ^/s = 314 GeV) require at least 
10 times less CPU time using partial parameterization. 
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6 Conclusions 



We have developed paramcterizations of electromagnetic showers for different materials 
and sampling geometries. Shower to shower fluctuations and correlations are taken into 
account consistently, as well as correlations between the longitudinal and radial shower 
development. Comparisons with data have shown that parameterized simulations are able 
to predict measured calorimeter signals with an acceptable precision. Using the methods 
described above, the energy resolution is reproduced at the level of d=0.5%. The energy 
deposited in longitudinal and lateral layers is predicted with a precision of typically ±1.5% 
for both, the means and the fluctuations. Using partial paramcterizations, the energy 
measured in electromagnetic (and hadronic) modules differs by an amount of 1.7% (9%), 
if the beam enters directly into a crack. The paramcterizations presented here provide a 
fast and precise algorithms for large scale Monte Carlo production of events for physics 
analysis. 



A Summary of formulae 

A.l Homogeneous Media 

A. 1.1 Average longitudinal profiles 



Thom = Iny- 0.858 

ahom = 0.21 + (0.492 + 2.38/Z)ln7/ 



A. 1.2 Fluctuated longitudinal profiles 



(hiThom) = ln(lny- 0.812) 

a{lnThom) = (-1.4 + 1.261ny)-i 

{Inahom) = In (0.81 + (0.458 + 2.26/Z) In y) 

a{\nahom) = (-0.58 + 0.86 In y)"^ 

p{lnThomMahom) = 0.705 - 0.023 In y 



A. 1.3 Average radial profiles 



Rc,hom{T) = Zi^Z2T 

RT,hom{T) = A;i{exp(/c3(r-/c2))+exp(A;4(r-/c2))} 



( \ P2- I 

Phom{r) = Pi exp <^ exp 



P3 V P3 
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with 

zi = 0.0251 + 0.00319 InE 

Z2 = 0.1162 + -0.000381Z 

A;i = 0.659 + -0.00309Z 

k'i = 0.645 

kz = -2.59 

ki = 0.3585 + 0.0421 In £; 

pi = 2.632 + -0.00094Z 

P2 = 0.401 + 0.00187Z 

P3 = 1.313 + -0.0686 In 

A. 1.4 Fluctuated radial profiles 

t exp((lna)) 



{t)i exp((lna)) - 1 

Tspot = Tftom(0.698 + 0.00212Z) 
aspot = a,,o^(0.639 + 0.00334Z) 

A. 2 Sampling Calorimeters 

A. 2.1 Material and geometry parameters 

Wi = — — [p = density) 

Pjdj 

i 



E 



^0,eff i Xo,i 

1 1 ^ ^ 21.2 MeV) 



RM,eff Eg . Xq 

i 



^0,1 



F. = 



e = 



da ~t~ dp 

1 



1 + 0.007(Zp - Za) 
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A. 2. 2 Average longitudinal profiles 

Tsam = Tw - - 0.53(1 - e) 

asam = Oihom " 0.444^^^ 

A. 2. 3 Fluctuated longitudinal profiles 

In (exp((ln7),„^)) - O.55F5 ^ - 0.69(1 - e)) 
(-2.5 + 1.25 In y)-^ 
In (exp((lnaw))-0.476F5i) 

(-0.82 + 0.79 In y)-^ 
0.784 - 0.023 In y 

A. 2. 4 Average radial profiles 

Rc,hom - 0.0203(1 - e) + O.O397F5 1 cxp(-T) 
RTMn^ - 0.14(1 - e) - O.495F5 1 cxp(-r) 
Phom + (1 - e)(0.348 - 0.642^5^ exp(-(r - 1)^)) 

A. 2. 5 Fluctuated radial profiles 

t exp((lnQ;)) 
(i)iexp((lnci;)) - 1 

T;,o^(0.813 + 0.0019Z) 
aw(0.844 + 0.0026Z) 



(lnr,„^) = 

o-(ln7;a„j) = 

(In asam) = 

(7(ln asam) = 



RCjSam 
RT,sam — 
Psam — 



Nspot = 
Tspot = 

OlSpot — 
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Figure Captions 





1 Depths of the shower maxima as a function of y (GEANT) 


2 Shape parameter a as a function of y (GEANT) 


3 Parameters of longitudinal profiles in homogeneous media (GEAN'i') . . 


4 Mean longitudinal profiles and their fluctuations in a lead glass calorimeter 


(GEANT and param. simulation) 




5 Parameters of individual longitudinal energy profiles in sampling calorime- 


ters (GEANT) 





6 Mean longitudinal profiles and their fiuctuations in the ill lead liquid argon 


calorimeter IFE (Hl-IFE) (GEANT and param. simulation) 





7 Energy containment vs. calorimeter depth in Hl-IFE (GEANT and param, 
simulation) 





8 Energy resolution vs. calorimeter depth in Hl-IFE (GEANT and param. 


simulation) 





9 Mean radial profiles in uranium at various calorimeter depth (GEANT). 


'I'he two component Ansatz is indicated for the depth 1 — 2Xq 






10 Shower depth and material dependence of radial shower parameters 





11 Tungsten, 40 GeV. Left figure: Decomposition of the two sources of radial 



shape fiuctuations. Right figure: Comparison of energy and spot density 


distributions 






12 Mean radial profiles and their fiuctuations in lead (40 GeV) and uranium 


(100 GeV) (GEAN'I' and param. simulation) 






13 'I'op: Mean integrated radial profiles in lead glass. Bottom: Correspond- 


ing shape fiuctuations: — (■ ■ ■) with (without) long./rad. correlations 



(GEAN'I' and param. simulation] 





14 'I'otal number of spots required to reproduce radial shape fiuctuations in 


sampling calorimeters] 




15 Mean radial profiles in Hl-IFE (GEANT and param. simulation) .... 


16 Radial shape fiuctuations in Hl-IFE (GEANT and param. simulation) 


17 Energy containment vs. calorimeter radius in Hl-IFE (GEANT and param. 




simulation) 






18 Energy resolution vs. calorimeter radius in Hl-IFE: — (■ ■ ■) with (without) 


long./rad. correlations (GEANT and param. simulation) 




19 Energy distributions of the cluster with the maximum number of channels 


{Emax) and of all other channels [Erest] for 30 GeV showers in three difi'erent 



modules, IFE, CB3, and FBI of the HI lead liquid argon calorimeter (data 



and H1FA5T) 

20 Mean longitudinal profiles and fiuctuations in Hl-IFE (data and HIFAST) 

21 Energy distributions in 4 subsequent layers of H1-CB3 for 30 GeV (data 
and HIFAST)! 
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22 Mean lateral profiles and fluctuations at ditt'erent depths in Hl-ll^'E for 10 


GeV (data and HIFAST). The coordinates J correspond to transverse pad 


numbers 




23 Mean lateral profiles and fluctuations in Hl-FBl (data and HIFAST). The 


coordinates 1 correspond to transverse pad numbers 




24 Normalized energies as measured across a crack between the modules HI- 


CB2 and H1-CB3 (data and HIFAST) 
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